Enhanced Molecular Orientation Induced by Molecular Anti-Alignment 
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We explore the role of laser induced anti-alignment in enhancing molecular orientation. A field- free 
enhanced orientation via anti-alignment scheme is presented, which combines a linearly polarized 
femtosecond laser pulse with a half-cycle pulse. The laser pulse induces transient anti-alignment in 
the plane orthogonal to the field polarization, while the half-cycle pulse leads to the orientation. 
We identify two qualitatively different enhancement mechanisms depending on the pulse order, and 
optimize their effects using classical and quantum models both at zero and non-zero temperature. 

PACS numbers: 33.80.-b, 02.30.Yy, 32.80.Lg 



I. INTRODUCTION 



The concept of molecular alignment refers to the an- 
gular localization of the symmetry axis of a molecule. 
The degree of alignment is measured by (cos^ 9) , where 
6 is the angle between the symmetry axis and a speci- 
fied space-fixed axis, and () denotes an average over the 
molecular ensemble. A molecule is said to be aligned if 
its symmetry axis lies along the space-fixed axis, whereas 
it is defined to be anti-aligned if its symmetry axis is per- 
pendicular to the space-fixed axis. In recent years, var- 
ious applications of molecular alignment have been pro- 
posed, including high harmonic generation 0], laser pulse 
compression , nanolithography 0] , control of photodis- 
sociation and photoionization 4] , and quantum informa- 
tion processing |5j|. These applications have stimulated 
the need to align molecules optimally under field-free con- 
ditions. An important milestone in the development of 
alignment methods is the use of linearly polarized, ultra- 
short laser pulses to create a rotational wave packet by 
an impulsive Raman mechanism. If the temporal pulse 
width of the laser is shorter than the rotational period of 
the molecule (i.e., if its bandwidth is greater than the ro- 
tational level spacing), the molecule undergoes a series of 
Raman excitations that produce a coherent superposition 
of rotational states For short pulses, peak field- free 
alignment along the electric vector of the laser field is 
achieved after termination of the laser pulse, at a time 
that depends on the pulse strength. As the wave packet 
evolves, the molecule loses its alignment, and even be- 
comes anti — aligned at some later time. The molecule 
also undergoes a series of field-free realignments at in- 
teger multiples of the revival time, Trev = l/(2i?c), where 
B = h/{A-KcIm) is the rotational constant, c is speed of 
light, and is the moment of inertia. In addition, a 
number of fractional rotational revivals occur at rational 
fractions of r^ev [El^' 

Molecular orientation refers to the case of molecules 
with a directional symmetry axis (i.e., when a molecule 
has a permanent dipole). Orientation is measured by 
(cos 9) , where 9 is the angle between the molecular dipole 



and the same spaced- fixed axis used to define alignment. 
Naturally, orienting a molecule requires breaking of the 
orienting field symmetry. Various symmetry-breaking 
methods have been proposed, including introduction of a 
weak DC electric (or magnetic) field in conjunction with a 
pulsed laser field , and coherent excitation with laser 
fields of frequencies u) and 2u; The most versatile 

method for orientating dipolar molecules utilizes asym- 
metric electromagnetic half-cycle pulses (HCPs) ■ 
Maximum field-free orientation is achieved some time af- 
ter application of the HOP. A single short HCP has a 
limited effect, however. As was shown in ref. its 
effect saturates with intensity. An HCP applied to a 
group of molecules that are initially oriented randomly 
in space contributes different angular velocities to indi- 
vidual molecules, so that molecules starting from obtuse 
angles {tt/2 < 9 < n) move too slowly to catch up with 
molecules starting from acute angles (0 < 6* < 7r/2). As a 
result, the kicked molecules do not all point in the same 
final direction at any given time, and the orientation pa- 
rameter has a maximum value of < cos 9 >~ 0.75. This 
effect is similar to non-perfect focusing caused by spher- 
ical aberration in geometrical optics. It was shown the- 
oretically 0, 0, 0| that the degree of field-free orien- 
tation and alignment can be enhanced by using trains of 
laser pulses, and enhanced alignment by a pair of pulses 
has been demonstrated experimentally |0 0, 0, |23| ■ 
(Fo r a recent review of field-free alignment, see Refs. 
j2ll [2^.') Creating a train of HCPs with sufficient 
strength to orient a molecule is experimentally difficult, 
however, requiring that alternative approaches be devel- 
oped. 

As we recently showed j2j] , anti-aligning molecules in a 
plane perpendicular to the HCP polarization can enhance 
the orienting effect of the HCP. In the current paper we 
investigate in detail the " orientation via anti-alignment" 
approach that combines an asymmetric half-cycle pulse 
with a symmetric femtosecond laser pulse. The latter 
induces anti — alignment, whereas the HCP orients the 
molecules. Depending on the pulse order, we identify two 
different mechanisms for enhanced field-free orientation, 
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and describe them qualitatively in section^ A rigorous 
formulation of the problem is given in 11111 A classical 
approach to the model at zero and finite temperatures 
is presented in section Hvl and a full quantum treatment 
is developed in section where we also demonstrate 
further enhancement via a three pulse scheme. Finally, 
we summarize our findings in section IVTl 



II. "ORIENTING AN ANTI-ALIGNED STATE" 
AND "CORRECTING THE ROTATIONAL 
VELOCITY ABERRATION" 

We identified two qualitatively different mechanisms 
for the orientation enhancement, which we term "orient- 
ing an anti-aligned state" and " correcting the rotational 
velocity aberration," depending on the temporal order of 
the pulses. 

The "orienting an anti-aligned state" mechanism is il- 
lustrated in Fig. ^ The top drawing depicts a ran- 
domly oriented collection of molecules. A short sym- 
metric laser pulse applied to the molecules pushes their 
symmetry axes toward the plane perpendicular to the de- 
sired orientation direction, preparing the molecules in an 
anti — aligned state angularly localized near 9 = ^/2. A 
schematic illustration of the anti — aligned state is given 
in the middle drawing. When a delayed strong asymmet- 
ric HCP is applied to such an ensemble, all the molecules 
gain nearly the same rotational velocity. Because all the 
molecules depart from the anti-aligned state with very 
similar initial angles, they reach the orientation direc- 
tion almost simultaneously at some later time, pointing 
in the same direction, as depicted in the bottom draw- 
ing. As is apparent from the bottom two drawings in Fig. 




FIG. 1: "Orienting an anti-aligned state" mechanism. (Top) 
The symmetric pulse is applied to an ensemble of randomly 
oriented molecules. (Middle) After a time ti, when the 
molecules are anti-aligned, an asymmetric orienting pulse is 
applied. (Bottom) After an additional time t2, the molecules 
are oriented. 

^ the more strongly the molecules are anti-aligned, the 
better they subsequently become oriented. 

The "correcting the rotational velocity aberration" 
mechanism is illustrated in Fig. |21 First, we apply the 
asymmetric HCP. The molecules gain angular velocity in 
the direction illustrated in Fig. |21 by the filled arrows. 



Then, shortly after the orienting HCP (or even simul- 
taneously with it), we apply the symmetric laser pulse. 
Such a pulse decelerates the rotation of molecular dipoles 
having acute angles with respect to the orientation direc- 
tion and accelerates those dipoles having obtuse angles. 
The directions of the angular velocity induced by the 
symmetric pulse is illustrated in Fig. [3 by the open ar- 
rows. This effect compensates for spherical aberration 
in the angular distribution of the rotational velocity dis- 
tribution and improves the overall orientation at a later 
time. 




FIG. 2: "Correcting the rotational velocity aberration" mech- 
anism. Two schematic molecules are plotted. The filled ar- 
rows denote the first HCP orienting direction, and the open 
arrows denote the second symmetric pulse directions. The 
molecule at S < 7r/2 is decelerated by the symmetric pulse, 
whereas the molecule at 9 > 7r/2 is accelerated by the sym- 
metric pulse. 

In the following sections we demonstrate for both 
mechanisms that significantly enhanced orientation may 
be achieved by a proper choice of the delay between the 
pulses and of their relative intensities. 

III. FORMULATION OF THE PROBLEM 

The Hamiltonian of a 3D driven rigid rotor (linear 
molecule) interacting with a linearly polarized field is 
given by 

H = ^ + vie,t), (1) 

where J is the angular momentum operator and 9 is the 
angle between the molecular axis and the polarization 
vector of the field. For a symmetric laser pulse interact- 
ing with the induced polarization, the interaction term, 
averaged over the fast optical oscillations, is given by 

V{9,t) = -y\t)[{a\i-ai_)cos^{9)+aA_], (2) 

where e{t) is the envelope of the laser field, and a|| and 
a±_ are the parallel and perpendicular components of the 
polarizability tensor, respectively. For a symmetric laser 
pulse the contribution from the permanent dipole aver- 
ages to zero. 

For an asymmetric HCP, the interaction with the 
dipole moment is given by 

Vi9,t) = -^e{t)cos{9), (3) 
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where /i is the permanent dipole moment and e(t) is the 
ampHtude of the HCP. In the present paper we assume 
that the duration of the laser pulse is much shorter than 
the contributing periods of the rotational wave packet, 
so that the excitation dynamics may be calculated in the 
impulsive limit. The impulse imparted to the rotator 
is characterized by a dimensionless kick strength (or ac- 
tion), P. For an asymmetric pulse, P is given by 

/oo 
e{t)dt, (4) 
-OO 

where the integration is performed over the unidirectional 
part of the HCP, and 

/oo 
e^{t)dt (5) 
-oo 

for a symmetric pulse. We start with the mechanism of 
"orienting an anti-aligned state" and consider a rotator 
excited first with a symmetric pulse of strength Ps at 
t = and then with an asymmetric pulse of strength Pa 
at t — ti. Henceforth the dimensionless time is measured 

in the units of Im/h =^ Trev/'^TT. 

IV. CLASSICAL TREATMENT 

Considerable physical insight may be derived from the 
(semi)-classical treatment of the problem, which is valid 
for Ps,Pa ^ 1- This is a natural approach to the process 
of enhanced orientation involving highly excited rota- 
tional states. In this section we formulate the orientation 
problem classically. A more general (although less intu- 
itive) quantum mechanical treatment of the same prob- 
lem is provided in sectional Classically, if a rotationless 
molecule is initially aligned at angle 9q, it will be found 
at the same angle just after the first symmetric kick (see 
Eqs.©,©) but with angular velocity — sin(26'o), so 
that at some later time it will have an angle 

eit) ^9o- Pstsm{29o). (6) 

When the orienting pulse of strength Pa is applied at 
time t ~ ti, the velocity increment is —PaSm[6{ti)], so 
that the angular velocity after the second pulse is Lu{da) — 
-Ps sin(26lo) - Pa sin[6'o - Psh sin(26'o)]. The angle 6 at 
time t2 after the second pulse is therefore 

e{ti+t2) = 9o- P,hsm{20o)~t2{PsSm{29a) 

+ PaSm[So-Pstism{2eo)]}. (7) 

A similar expression may be derived for the inverse 
order of pulses, utilized in the second mechanism of en- 
hanced orientation. The alignment at time t = ti and 
orientation at time t = ti + 12 is calculated by averaging 
cos*^ d{t) over all values of Oq, 

1 

(cos*-' 0{t)) = - / cos'' 9{t) sin eodOo, (8) 
2 Jo 



where k = 1 and 2 for orientation and alignment, respec- 
tively. 

Our analysis shows that strong transient anti- 
alignment may be achieved via two related methods. The 
first one is of a classical nature. It operates on a short 
time-scale (t << Trev) and requires negative values of 
the kick strength P,- Pulses with negative P, produce 
anti-alignment by pushing molecules into the equatorial 
plane {9 = tt/2). Fig. [SJl shows the expectation value 
of (cos^ 9) calculated both classically and quantum me- 
chanically (as will be explained in section O for a pulse 
with Pg = —10. Both treatments predict a deep min- 
imum ((cos^ 6'),„i„ = 0.077) shortly after the pulse at 
trn ~ O.S/jPsl. There are various ways of achieving a 
negative kick strength of the symmetric pulse. First, 
some alkali halides (such as LiF) have a negative polar- 
izability anisotropy (a|| — a i ) |24| . Second, the interac- 
tion of molecules with a circularly polarized light pulse 
propagating in the direction of the desired orientation is 
proportional to P, sin^ 9 — Pg ~ Ps cos^ 9. This result is 
formally equivalent to the interaction with a linearly po- 
larized pulse (our model) having negative kick strength. 

A third method of achieving ant i- alignment uses laser 
pulses with positive Pg, which cause a substantial 
alignment on a short (classical) time scale after the kick 
(see Fig. |31d). However, if it were possible to invert the 
dynamics and travel backwards in time, one would ob- 
serve a strong anti-alignment before the kick. Remark- 
ably, the effect of quantum revivals provides such an 
option. Indeed, the rotational wave function satisfies the 
following equation 

^{9,Trev-t)=-i'{9,-t), (9) 

and, for sufficiently strong pulses, quantum dynamics 
just before one full revival cycle is equivalent to classical 
dynamics analytically continued to negative times. As 
a result, considerable ant i- alignment is observed in this 
time domain (see Fig. ^jp). 

Next we would like to optimize the orientation fac- 
tor \{cos{9)) {Pg , Pa, ti, t2) \ using the classical model (Eqs. 
ITI) and ©) with an anti- aligning pre-pulse {Pg < 0). In 
this procedure, the model is formally extended to nega- 
tive times to cover effects in the full revival time-domain, 
as explained above. At zero initial temperature, the op- 
timal solution depends only on the ratio Pa/ Ps and the 
products Pgti and Pat2- A contour plot of the orienta- 
tion factor as a function of |Ps|ti and \Pg\t2 (for a fixed 
ratio Pa/\Ps\ — 3) is given in Fig. ^ where two almost 
equivalent extremum points are observed. 

Figure |S1 displays the highest (optimized with respect 
to the delay time) post-pulse orientation parameter, 
(cos 9) (solid line), and the optimal delay between pulses, 
tiopt (dashed line), as functions of Pa/|Ps|. There are two 
optimal solutions of almost the same efficiency. The first 
one provides maximal orientation in the direction of the 
HCP shortly after the second pulse (upper panel in Fig. 
|3J). The second one (lower panel in Fig. [SJ delivers en- 
hanced orientation in the opposite direction in the full 
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time 



FIG. 3: Alignment parameter vs time after excitation by a 
laser pulse with (a) Pa = —10 and (b) Ps = 10. Solid curves 
are the quantum results. Dashed curves are calculated clas- 
sically (a) for positive time and (b) for negative time (shifted 
by Trev)- Strong anti-alignment is seen in both cases. 




FIG. 4; Classical orientation factor as a function of |Ps|ti 
(horizontal axis) and |Ps|t2 (vertical axis) for the first mecha- 
nism (the laser pulse precedes the HCP). Pa/\Ps \ ~ 3. There 
are two almost equivalent extreme points. 



revival domain after the second pulse. In both cases, 
an impressive value of [(cos^)! « 0.95 is achieved for 
rather modest pre-pulses {\Ps\ ^ O.IPq) (as compared 
to the limit of {cos6)max ~ 0.75 for a single HCP). In 
this regime, the optimal delay between pulses, i.e. opti- 
mal ti, asymptotically approaches the time of the best 
anti-ahgnment, tm ~ 0.8/|Ps|. 

It is easy to show that the same degree of orientation 
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FIG. 5: Classically optimized orientation factor (solid curves) 
and delay between pulses (dashed curves). The laser pulse is 
fired before the half-cycle pulse. 

may be achieved by combining an HCP with an aligning 
laser pulse {Pg > 0). The apparent symmetry relations, 

{cos 9) (Ps, Pa, -tl,-t2) = {C0se)i-Ps,-Pa,h,t2) 

(10) 

and 

{cOs9){Ps,Pa,h,t2) ^ -{cOs9){Ps,-Pa,ti,t2), (11) 

reduce this problem to the already-studied case of the 
anti-aligning pulse. 

We used the same approach to analyze the second 
mechanism of enhanced orientation mentioned in the In- 
troduction. In the simplest, non-optimized version, the 
orienting and anti-aligning pulses are applied simultane- 
ously {ti = 0). Direct numerical maximization of the 
expression in Eq. |S1 shows that {cos9)max — 0.89 when 
Pa/\Ps\ ~ 2.34 and \Ps\t2 ~ 0.78. When the hybrid 
pulse is composed of an orienting component and an 
aligning one {P^ > 0), the symmetry relations (|10|l and 
(|ll|l predict the same orientation, but in the opposite 
direction {{cos9) = —0.89) just before one full revival 
cycle {Pst2 ~ —0.78). This effect was reported in a re- 
cent paper |25l | as a result of direct numerical simulation 
of the quantum rotational dynamics of molecules excited 
by a single hybrid pulse. More efficient results may be 
obtained when the HCP precedes the anti-aligning laser 
pulse. In this case, our classical analysis reveals a sin- 
gle dominating optimal solution. A contour plot of the 
orientation parameter is given in Fig. The optimal 
dominant solution is presented in Fig. [T] The maximum 
orientation {cofi9)„iax ~ 0.96 is reached at Pa/\Ps \ ~ 1.6, 
and the optimal delay is tm ~ 0.36/|Ps|. 

As the next step, we calculate the orientation factor 
at a nonzero temperature. We shall derive the orienta- 
tion factor for the " orientation of an anti-aligned state" 
mechanism only; the orientation factor for the second 
mechanism may be similarly derived. 
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FIG. 6: Classical orientation factor as a function of [Pslti 
(horizontal axis) and |Ps|i2 (vertical axis) for the second orien- 
tation mechanism (HCP precedes the laser pulse). Pa/\Pa \ ~ 
1.6. There is a single dominant maximum. 




FIG. 7: Classically optimized orientation factor (solid curve) 
and delay between pulses (dashed curve). The laser pulse is 
fired after the half-cycle pulse. 



The Lagrangian for the 3D free rotor is given by 



L = (9^ + 6^ sin^ ( 



and the Euler-Lagrange equation for 9 is 
d dL dL 



dt d9 d9 



0. 



(12) 



(13) 



The canonical momenta are 
dL 



P. 



I(j> sin^ 9, 



(14) 



(which is a constant of motion as </> is a cyclic coordinate), 
and 



(15) 



From Eqs. l(T^ . ifT^ . and ((T5|) the following equation of 
motion follows: 



-P| cos 9 



P 



sm 



(16) 



In what follows, it is convenient to measure the canon- 
ical momenta in the units of pth = Im^th with ujth = 
y/ksT / I„i, where T is the temperature and /cs is the 
Boltzmann constant. By setting P'^ = P^/pth,Pg = 
Pe/pth, and t' = cotht, one writes Eq. as 



9^ P: 



:I2 



COS 9 
sin^ 9 ' 



(17) 



The solution to this equation is 

1,. PL 



cos 9{t') 



LU 
P' 



cos(6'(0) - ujt') 



)cos(6'(0) H-cji'), (18) 



where 



P'l 



a/2 



sin^ 9{Q) ' 



(19) 



Here P'g is a constant initial canonical momentum (see 
Eq. iQSJ). As a result of a kick induced by a symmetric 
laser pulse applied at t = 0, Pg changes its value to 



P^=P^(O)-P;sin(20(O)). 



(20) 



The initial thermal distribution function for the 
molecules has the following form in the chosen dimcn- 
sionless variables: 



/(0, 0, p^, p;) = ^ cxp[- i(pf + ^)]- 



sm 



(21) 



The thermally averaged alignment factor reads as 

pTT nllZ nOG poo 

{cos' 9){t') = / d9{0) dm rfP^(O) / dP'iO) 











X cos' 9{t') f (9(0), m,PeiO),P^m, (22) 

where / is given by Eq. (|21|l . The factor cos 9{t') is 
provided by Eq. in which P^ = P^(0), and P^ is 

given by Eq. (|20|l . An HCP is applied with a delay t[ 
after the first symmetric laser pulse. Using the equations 
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P'g=9 and d{cose)/dt 
second pulse: 



sin 99 we obtain Pn after the 



d(cos0)M_ .^^^ 



sm( 



(23) 



Here sin 9 and cos 9 are obtained ft-om Eq. H18(l at t' = t[ . 
The new value of Pg and the angle 9{t[) are substi- 
tuted again into Eq. (|18|l to propagate the rotation angle 
further. Thus, the time-dependent orientation factor is 
given by 



{cos 9) {t[+t'^: 



X cos9{t' = t[ +t'2) 

X f{9iO),m,PoiO),P^m- (24) 

From here on, we return to the previous dimensionless 
units used in this paper. In particular, the dimensionless 
kick strengths are replaced by 



^^(O) / d(t>{0) 
Jq 

dpm / dP'M 



p' 



Pn 



(25) 



Figure IHl displays a comparison between quantum (see 
the next section) and classical calculations for KCl 
molecule kicked by two coinciding pulses {ti — 0) at 
T = 5K. One may observe that the classical and quan- 
tum results agree well for strong pulses (P 3> 1) and 
short propagation time. 
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FIG. 8: Orientation parameter vs. t2 {ti = 0) for T = 5K, 
Pi — —10, and P2 = 50. Solid and dashed lines are quantum 
and classical calculations, respectively. 



V. QUANTUM TREATMENT 

Consider a quantum rotator being initially in the 
ground state and kicked by a short symmetric laser pulse 



at t = 0. In the impulsive approximation, its wave func- 
tion acquires a phase factor as a result of the kick: 



-^{9, t = +0) = exp[iPs cos^ 9]Y^{9). 



(26) 



By expanding this expression as a sum of spherical har- 
monics, one finds the wave function at later time t: 



t) = ^ cj exphz J(2 J + l)t]Y2°ji9). (27) 
The coefficients cj are given by 0| 

X iFi[J +1/2, 2J + 3/2, iP^l (28) 

where is the confluent hypergeometric function. At 
t = ti the rotator is kicked by the orienting pulse, ac- 
quiring an additional phase factor 

-^(9,11+0) = exp[iPaCos9]'^i9,ti - 0). (29) 

We use the well-known expression 

00 

exp(zP, COS0) = *^V47r(2J+l)j,(^a)r/(0), (30) 
,7=0 

where jj{Pa) is a spherical Bessel function, and again ex- 
pand the wave function in a series of spherical harmonics 



vi/((?,ti + o)-^diyA^), 



(31) 



(=0 



where 



di = -4f= f] E*V(2J+ l)ji(-Pa)Q' exp[-tl'{2l' + l)h] 

I' =01=0 



X y/{2j + 1){AI' + 1)121 + 1) 



C{j,2l',l\0,0,0y 
21 + 1 



(32) 



Here C{j, 21' ,l\0, 0,0) is a Clebsch-Gordan coefficient. 
This new wave function is allowed to propagate freely 
until t = ti + t2, at which point the orientation and 
alignment parameters are calculated by 

(cos'^e) = {9, t)\ cos'' 9\^ {9, t)),k^ 1,2. (33) 

We performed a fully quantum-mechanical analysis of 
the "orientation via anti-alignment" and "correcting the 
rotational velocity aberration" mechanisms, using the 
above-described methodology. Figure IHl demonstrates 
the optimized values of the anti-aligning pulse strength, 
\Ps\, the optimum delay between pulses, tm, and the 
maximal orientation parameter, {cos 9)max, as function 
of Pa- Very good agreement between quantum and 
classical results is observed even for moderate anti- 
aligning and orienting pulses {Ps,Pa ~ 3). Remarkably, 
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a significantly enhanced orientation may be achieved 
with field strengths available currently in the laboratory. 
Considering a KCl molecule in the ground state (having 
a revival time trev ~ 128ps, a dipole moment fj, ~ 10. 3D, 
and a polarization anisotropy — a±) « 3.1A , data 
taken from [l^l), one expects Pa ~ 10 for an HCP with 
the amplitude of SbkV/cm and a duration of about 2ps 
(1/e half width). According to Fig. El the orientation 
factor (cos 6') w 0.95 may be observed if the HCP is 
followed by a delayed anti-aligning pulse of 2ps duration 
and 5 x IQ^^W/cm? peak intensity. Negative times 
in Fig. |51 correspond, as was explained above, to the 
dynamics in the full revival domain. 




a 




2 4 G S 10 I 2 4 B e 10 



Pa Pa 



FIG. 9; Quantum mechanical optimized results at zero tem- 
perature. The left column (a,b) corresponds to the mech- 
anism of "orienting an anti-aligned state" (laser pulse fired 
before the HCP). The right column corresponds to the mech- 
anism of "correcting the rotational velocity aberration" (in- 
verse order of pulses). Upper panels (a,c) display the optimal 
strength of the anti-aligning pulse (solid lines) and delay be- 
tween pulses (dashed lines) as a function of the HCP strength. 
Lower panels (b,d) present the highest value of the post-pulse 
orientation factor vs the strength of the HCP. 

We also analyzed quantum mechanically the orienta- 
tion factor after two-pulse excitation for an ensemble hav- 
ing a non-zero temperature. The details of calculations 
can be found in the Appendix to this paper. The opti- 
mized solution for the orientation of an anti-aligned state 
is plotted in Fig. ^|as a function of Pq. A comparison 
between the classical and quantum thermal results was 
given in Fig. |S1 One may observe in Fig. ^| that the 
temperature reduces the orientation effect. However, the 
latter is still considerably better compared to the single 
HCP excitation. The maximal orientation factor for KCl 
excited with a single HCP of strength Pa — 10 is 0.59 at 
T = 5K. The combined action of the HCP and laser 
pulse provides a much higher value of 0.69 at the same 
thermal conditions. To achieve even better orientation, 
one may try to add additional laser pulses to improve the 



FIG. 10: Quantum mechanical optimized results for KCl at 
T = 5K. The graphs correspond to the mechanism of orient- 
ing an anti-aligned state (laser pulse fired before the HCP). 
a) Anti-aligning pulse strength and delay between pulses vs. 
HCP strength are given by solid and dashed lines, respec- 
tively, b) Orientation parameter vs. HCP strength. 



anti-alignment in the system. The simplest scheme of this 
kind consists of two symmetric laser pulses separated by a 
delay ti , followed by an HCP after a second delay, t2 ■ The 
optimal field-free orientation occurs at some time after 
the HCP. A full quantum optimization of this scheme at 
finite temperature presents a challenging numerical task, 
and goes beyond the framework of the present paper. The 
results of optimization at zero temperature are shown in 
Fig. ^2 In particular, using the experimentally feasible 
HCP discussed above (i.e. Pa — 10 at 85^:^/ cm) we ob- 
tain the optimized orientation factor (cos(0)) = 0.97 for 
Psi = -1.42 and Ps2 = -3.95. This means a 25% reduc- 
tion in the width of the angular distribution compared 
to the two-pulse orientation scheme. 



VI. CONCLUSIONS 

We presented an experimentally feasible method for 
enhanced orientation of linear dipolar molecules using a 
half-cycle pulse combined with a delayed laser pulse in- 
ducing molecular anti- alignment. Two qualitatively dif- 
ferent enhancement mechanisms were identified depend- 
ing on the pulse order, and their effects were optimized 
with the help of quasi-classical as well as fully quan- 
tum models. The transparent physics behind the in- 
terplay between anti- alignment, orientation, and quan- 
tum rotational revivals provides a solid basis for the 
future design of more sophisticated and efhcient solu- 
tions. In particular, we demonstrated that enforced anti- 
alignment by a pair of symmetric laser pulses prior to 
applying an HCP improves the orientation even more. 
These results may be generalized to trains of multiple 
laser pulses, similar to forced multi-pulse alignment tech- 
niques m [11 El El m El Ellg . In the same way, 
the "spherical aberrations" of a single HCP may be bet- 
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(a). 




In Eq. (PI, 

^rw = (cwico«(^)ic°w) (36) 

is the the contribution to the orientation parameter from 
initial state \lQ,mo). 

The wave function \ipj^"{t)), is calculated as follows. 
After the first symmetric laser pulse, the wave function 
becomes 

IC°(0^)) =exp[zP,cos2(0)]|/o,mo). (37) 

Next, we expand {"ipj^^iO^)) in the spherical harmonic 
basis. 



FIG. 11: Quantum mechanically optimized results for the 
three-pulse mechanism, where the first two pulses are sym- 
metric, and the third pulse is an asymmetric HCP. a) First 
pulse strength Pi and the delay between first and second 
pulses ti vs. HCP strength are given by solid and dashed 
lines, respectively, b) Second pulse strength P2 and the de- 
lay between second and third pulses vs. HCP strength are 
given by solid and dashed lines, respectively, c) Orientation 
parameter vs. HCP strength. 



ter corrected by multiple series of delayed laser pulses 
as well. Finally, a hybrid pair of delayed pulses may be 
followed by a series of well-timed symmetric laser pulses 
designed to preserve the achieved orientation over an ex- 
tended time period |26L . 
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C(0^) = E"ri^™o), (38) 

lo=0 

where the quantum number mg is preserved. The coeffi- 
cient a™" may be derived by expanding exp[iPs cos^(0)] 
in a sum of spherical harmonics using Eas. (|26l) . 127|l . and 
(|28|l and applying the rules of angular momentum alge- 
bra. However, we employ here another mathematically 
equivalent approach that is numerically advantageous for 
large values of P,- For this, we expand exp[iPs cos^{9)] 
in a series of Legendre polynomials, using the formula 
similar to Eq. (|30() . which gives (for negative Pg): 

expltPs cos^e)] = exp[-z^] z3^(2J + l)jj(^) 



,7=0 



X Pj(cos(26l)). 
Next, we express P,7(cos(20)) as 



2 
(39) 
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support of this research, and RJG acknowledges the Mo- 
torola Corporation and the US Department of Energy for 
its support. 



VIII. APPENDIX 

We present here the finite-temperature results for the 
"orientation of an anti-aligned state" mechanism only. 
The orientation factor for the second orientation mecha- 
nism may be derived similarly. 

At T > 0, the orientation factor should be thermally 
averaged over all initial states, so that 

(cos0) = ^P(?o) E ^rW' (34) 

Iq=0 mo = -lo 



where 



p„ , 1 r ^0(^0 + l) 



(35) 



th 



is the thermal distribution at temperature T, Q is the 
rotational partition function, and ath = {ksT /2Bhc)^/^ . 



Pj(cos(20)) = VrfL,jPL (cos 



L=0 



(40) 



where the coefficients d^.j are generated by an efficient 
recurrent procedure described in Appendix A of Ref . |3 . 
Using the well-known relation between Legendre polyno- 
mials and spherical harmonics, one finds al"" to be 



a. 



00 00 

"'sh 



expH^]2^2.^^.^^ (2-^+1) 



J=0 L=0 



. |P.| 2l„ + l 

X C{L,loJ\0,0,0)CiL,lo,l\0,mo,mo), (41) 

where C{L,lo,l\0,0,0),C{L,lo,l\0,mQ,mo) are Clebsch- 
Gordan coefficients. Free evolution after the first pulse 
gives 

00 

C(ir) = E"rea;ph^K^ + l)ii]l^,"^o). (42) 

lo=0 

Next, we apply the second (asymmetric) pulse at time ti, 



C(i^)=exp[zPaCos(0)]V,7(ir), 



(43) 
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and again re-expand the new wave function, 

oo 

V'r(^^)=E"ri^''"o), (44) 



in order to find the spherical harmonic coefficients: 



oo oo 



J=0 1=0 



X ^J^^C{J,lJ\0,0,0)C{J,lJ\Q,mo,mo). (45) 



Finally, this new wave function propagates freely for a 
time interval t2, giving 



V',7(i = h+t2) = ^ap°exp[--r(r+l)t2]|r,mo). (46) 

1=0 



This last expression is substituted into Eq. (|36|1 in order 
to calculate the orientation factor given by Eq. 1)34(1 . 
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